Strongly Interacting Fermions and Phases of the Casimir Effect 
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With the intent of exploring how the interplay between boundary effects and chiral symmetry 
breaking may alter the thermodynamical behavior of a system of strongly interacting fermions, 
we study the Casimir effect for the setup of two parallel layers using a four-fermion effective field 
theory at zero density. This system reveals a number of interesting features. While for infinitely 
large separation (no boundaries), chiral symmetry is broken or restored via a second order phase 
transition, in the opposite case of small (and, in general, finite) separation the transition becomes 
first order, rendering effects of finite size, for the present setup, similar to those of a chemical 
potential. Appropriately moving on the separation-temperature plane, it is possible to generate 
a peculiar behavior in the temperature dependence of the thermodynamic potential and of the 
condensate, compensating thermal with geometrical variations. A behavior similar to what we find 
here has been predicted to occur in bilayer graphene. Chiral symmetry breaking induces different 
phases (massless and massive) in the Casimir force separated by critical lines. 



It is well known that the presence of boundaries in 
empty space causes deformations of the quantum vacuum 
producing a macroscopic force. Casimir was the first to 
notice this for the case of two parallel, perfectly con- 
ducting plates and for the electromagnetic field, show- 
ing that the force, F c , is attractive. In natural units, 
F c = — 7r 2 S'a~ 4 /240, with a being the plate separation 
and S their surface area 

Since its original discovery, the study of the Casimir ef- 
fect has branched out in several directions and motivated 
experimental and theoretical research in an attempt to 
uncover its fundamental importance @, Q • A somewhat 
special attention has been directed towards understand- 
ing how the magnitude and, especially, the sign of the 
force depend on the properties of the vacuum, the geom- 
etry of the boundaries, or the external conditions, and the 
possible technological implications of this are by now well 
appreciated. Regarding this point, the fermion Casimir 
effect has triggered some initial curiosity due to the dif- 
ferent statistics obeyed by fermions as opposed to bosons, 
but it was soon understood that this difference does not 
necessarily lead to a change in the sign of the force. The 
application of the fermion Casimir effect that is most cel- 
ebrated is, perhaps, the bag model of hadrons (see Ref. [H 
for a review), but the scrutiny of fermion quantum vac- 
uum fluctuations has included the study of geometrical, 
thermal, dimensional effects, and others (see, for exam- 
ple, Refs. [H-UH or Ref. for a longer list of references). 

All previous work in relation to the fermion Casimir 
effect has focused, as far as we are aware, on the case of 
free fields. However, the problem may become more in- 
teresting if interactions are switched on. In fact, while in 
the case of free fields the dynamics is somewhat trivial, 
the inclusion of interactions opens up onto a richer spec- 
trum of possibilities. This problem is worthy of attention 
from different perspectives. First of all, the interplay be- 



tween chiral symmetry, the geometrical and topological 
properties of the spacetime, and quantum vacuum fluctu- 
ations is a concrete example allowing one to explore how 
the phenomena of symmetry breaking and the Casimir 
effect are connected (see, for example, Refs. 16- 2fjj|). 
Secondly, strongly interacting fermions have a well rec- 
ognized importance in describing several nonrclativistic 
condensed matter systems routinely studied in the labo- 
ratory. These include superconductors, conductive poly- 
mers, and several carbon-based materials (see Ref. [21| 
for some examples). Also, since chiral symmetry resides 
in the quark sector, strongly interacting fermion field the- 
ories are central tools to describe the phase diagram of 
QCD (see Refs. [22 , 23 for review). In this context, quan- 
tum vacuum energy effects have a direct relevance for any 
analysis of the phase structure aimin g at including effects 
of finite size (see, for instance, Ref. [24j). 

Let us begin our discussion summarizing a few points 
regarding the Casimir effect for free fermions. In the 
following, we will consider the case of parallel plates. 
Natural units will be used. The free fermion dynam- 
ics is described by the Dirac equation with the fields 
forced to obey some boundary conditions at the plates 
that we assume to be located at z = and z = a. 
Imposing, for simplicity, bag boundary conditions, ex- 
pressed as (1 + ij z ) i^\ z=0 a ~ 0, leads to an implicit 
constraint for the momenta in the z direction: 3?(fc z ) := 
msui(k z a) + k z cos(k z a) = 0, where m is the mass of the 
fermions [2|. The regularized Casimir energy, after inte- 
gration over the unconstrained directions, can be written 
as 
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where s is a regulator, fi is a renormalization scale, and 
the sum is over the roots of $(fc 2 ) = 0. In the massless 
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case, m = 0, one finds S = -77r 2 a" 3 / 2880 - When the 
fcrmion mass is nonzero, the above expression ([1]) for the 
Casimir energy can be recast in the following form, 
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where we have introduced the dimensionless quantity 
£ = ?na. The dots represent terms that do not con- 
tribute to the force. For £ = 0, expression ([2|) repro- 
duces the massless result given above. When £ ^ 0, the 
integral in ^ can be evaluated by expanding the inte- 
grand in the region of interest. For £ ^> 1, the force is 
exponentially suppressed. Basically, £ works as a mod- 
ulating parameter for the Casimir energy or force. For 
free fields, however, modulations caused by a change in 
the mass cannot occur dynamically (neither at zero nor 
at finite temperature or density) since, in the absence 
of interactions, symmetry breaking does not occur, and 
the mass is set by the chiral symmetry at the level of 
the Lagrangian [6j. On the other hand, when fermions 
are strongly interacting, chiral symmetry breaking occurs 
dynamically. This generates a mass for the fermions and 
is expected to induce a phase transition in the Casimir 
force. In a Casimir-like setup, the problem becomes par- 
ticularly amusing, since chiral symmetry breaking can be 
triggered not only by thermodynamical effects, but also 
by changes in the geometry or topology of the system. 

Such a mechanism is reminiscent of the mass genera- 
tion phenomena for self-interacting scalars in topolog- 
ically nontrivial spacetimes discussed, for instance, in 
Refs. [TMl. There, it was a combination of scalar self- 
interactions and nontriviality of the geometry responsible 
for the spontaneous symmetry breaking. Here, the mech- 
anism is driven by the breakdown of chiral symmetry due 
to thermal and geometrical effects, and it is controlled by 
the appearance of a condensate. 

To illustrate this idea with a concrete computation, we 
will consider the following theory: 



J2f = ^d^j + X (^) 



2N 



(3) 



At tree level, fermions are massless and the action is in- 
variant under discrete chiral symmetry, ip — > 7 s -0. N 
represents the number of fermion degrees of freedom and 
g is the coupling constant. In the following we will use the 
path integral approach and introduce finite temperature 
effects by means of the Matsubara formalism, t —> it with 
r e S 1 of period j3 = 2tt/T, where T is the temperature 
and antiperiodicity conditions imposed on the fermions. 
The presence of the condensate, a = — g (iJip) / 'N , can 
be made explicit in the partition function, Z, using the 
Hubbard-Stratonovich transformation, 



Z = I S>[4>, ip, cr] expf i I drdvD^ e Q 



where Jzf e g- = tfji^^d^^ — j^<t 2 — atpip. In the largc-iV 
approximation, an expansion of the path integral gives, 
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The above functional determinant has to be computed 
consistently with the boundary conditions imposed at 
the plates. Approximate analytic expressions can be ob- 
tained in specific regimes of temperature and separation 
251 ] . Here, we prefer to use a fully numerical approach 
that has the advantage of being more expedient and read- 
ily applicable to all parameter ranges. This adopts the 
method described, for example, in Refs. (26l - [28| to con- 
struct an appropriate contour integral representation for 
the functional determinant in @, that we evaluate nu- 
merically after a convenient contour deformation. Di- 
vergences are dealt with by means of zeta function reg- 
ularization, and finite temperature summations are car- 
ried out numerically. The theory ([3]) in four dimensions 
is nonrenormalizable and requires the introduction of a 
cutoff scale. Here, we will fix the coupling constant and 
the renormalization scale to achieve a broken symmetry 
phase at zero temperature in the absence of boundaries. 
Since the effective potential rescales as ft — > X D+1 fl un- 
der mass redefinition, mass — > A mass, (see Sec. Ill of 
Ref. HH), we may fix, without loosing generality, the 
value of the renormalization scale (numerical value is set 
to fi = 100) and express all quantities accordingly. 




Figure 1. Temperature dependence of the effective potential 
for fixed separation. The left-hand plot refers to the case of 
infinitely large separation for which the phase transition is 
second order. The curves refer (bottom to top) to the values 
of T/T crit = 0.00, 0.60, 0.80, 1.00, 1.17. The right-hand 
plot refers to the case of small separation a x T cr u = 0.6 in 
which case boundary effects are non-negligible and the phase 
transition becomes first order. The curves refer (bottom to 
top) to T/T crlt = 0.00, 0.76, f .00, 1.13, 1.29. 

Our results arc summarized in Figs. [TJH Figure [1] 
shows how the effective potential Q (normalized by sub- 
tracting its value Qq at a — 0) changes with temperature 
when the separation a is fixed. The right- (lcft-)hand 
panel refers to small (large) values of the separation a, 
and chiral symmetry breaking occurs via a first (second) 
order phase transition. Figure [2] illustrates the oppo- 
site situation, that is, the dependence of the effective 
potential on the separation a for fixed temperature. The 
left-hand panel of Fig. [2] refers to the case of very small 
temperature, while in the right-hand panel we have set T 
close to its critical value in absence of boundaries. From 



a a 
Figure 2. Separation dependence of the effective potential 
at fixed temperature. The left-hand plot refers to the case 
of small temperature (a CT u x T = 0.003) and the phase tran- 
sition is first order. The curves (bottom to top) refer to the 
following values: a/a crlt = 100, 2.79, 1.39, 1.00, 0.83. In the 
right-hand plot we have set the temperature close to its crit- 
ical value in the absence of boundaries and the curves refer 
(bottom to top) to a x T crit = 247.00, 7.41, 3.70, 2.47, 1.85. 

Figs. [1] and [21 it is clear that a decrease of a tends to in- 
crease the effective temperature of the system. However, 
when the separation is small and both geometrical and 
thermal effects arc important, the phase transition is a 
first order one, in contrast with a second order transition 
expected when changes are purely thermodynamical and 
boundary effects negligible. In this sense, effects of finite 
size for the present geometry seem to be more similar to 
those of a chemical potential than of temperature. 

The nature of the change in the order of the transition 
can be better understood using the Ginzburg-Landau ex- 
pansion for the thermodynamic potential, 

ft - ft = co(a, T)a 2 + c x (a, T)a 3 + c 2 (a, 1> 4 + ■ • • , 

where the coefficients Cj(a, T) are dimensionful functions 
of the temperature and separation and depend on the ge- 
ometry, topology, and the boundary conditions. In the 
absence of boundaries, i.e., for a —> oo, chiral symme- 
try, a O — a, prohibits the appearance of odd powers 
in the above expansion and guarantees the coefficients 
associated with such terms to vanish, leading to a sec- 
ond order phase transition. In the present case, how- 
ever, dimensional inspection indicates the presence of 
terms proportional to a a , and explicit computation 
of the associated coefficient shows that they do not van- 
ish (see Ref. [2j|). This causes the transition to change 
to first order. A deeper connection can be drawn us- 
ing the Schwinger-De Witt expansion that relates the 
coefficients c,(a,T) to the heat- kernel coefficients $i/2 
associated with the operator in Q, a(a,T) oc Q i j 2 (see 
Ref. [HI)- In this way, it is possible to see that odd pow- 
ers of the condensate in the Ginzburg-Landau expansion 
are accompanied by half-integer (boundary) heat-kernel 
coefficients. These (see Ref. [26( for explicit expressions 
of 9i/2 and Ref. [3(| for a clear discussion in relation to 
the Casimir effect) are related to the boundary geome- 
try, topology, and also to the boundary conditions that, 
in the present situation, break chiral symmetry. 

The above results suggest that simultaneous changes 
in temperature and separation may have interesting ef- 
fects. In fact, if a decrease in temperature typically 
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Figure 3. The figure shows how the effective potential changes 
when both temperature and separation decrease linearly ac- 
cording to T(S) = u — vS and a(5) = q — pS where u, 
v, q, p are constants and 8 is varied. In the left- (right- 
)hand panel we have set u = q — 2 and v — p — 1 
(u = p = 1 and 1/v = q = 30). The curves 1-6 in 
the left- (right-)hand panel correspond, respectively, to the 
values of 5 = -0.5, 0.095, 1.15, 1.55, 1.638, 1.7 (S = 
-60., -44., 15., 29.4, 29.64, 29.72). 
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Figure 4. aT phase diagram. The thick dashed (red) curve 
represents the critical line (the blue dots superposed are calcu- 
lated numerically) . The dots connected by the straight dotted 
brown (dot-dashed green) line refer to the values of a and T 
used in the left- (right-)hand panel in Fig. 

pushes down the potential and tends to bring the system 
into a broken phase with a nonvanishing condensate, a 
decrease in distance has the opposite effect. Therefore, 
by moving appropriately on the (a-T) plane, it is pos- 
sible to compensate thermal with geometrical changes. 
This may result in a peculiar phase diagram with the 
system moving towards a broken chiral symmetry phase 
despite an increase in temperature. An example of this 
sort is shown in Fig. [3] To illustrate this point, let us 
first look at the left-hand panel of Fig. [3l The initial 
configuration (curve 1) is obtained by setting T large 
and a small. In this case, the effective temperature is 
high and the system is in a symmetric phase. We then 
start decreasing both temperature and separation. (For 
illustration, we vary a and T linearly, according to the 
relations T(S) — u — vS and a(<5) = q — p5, where we 
set u = q = 2, v = p = 1 and change 5). The initial 
configuration (curve 1) refers to S = —0.5. Increasing S 
(i.e. lowering temperature and separation), the effective 
potential moves downwards and we encounter a second 
order phase transition for S ~ 0.095 (curve 2). A further 
increase in 5 keeps pushing the potential down and the 
condensate up until a reaches a maximum (for S ~ 1.15, 
curve 3). In this process, the breakdown of chiral sym- 
metry is driven by the temperature, assisted by effects 
of finite size that are responsible for the change in the 
order of the transition. From this configuration, a fur- 
ther increase in S changes the tendency of the potential 
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that, despite the decrease in temperature, starts to move 
upwards (and in a that descends from the maximum). 
When T and a reach a critical value (for 8 ~ 1.638, curve 
5), another second order transition occurs restoring chi- 
ral symmetry. In this region, it is the separation that 
dominates over the temperature. Starting from a config- 
uration where the initial separation is large, it is possible 
to reduce the effects of the boundaries adjusting the order 
of the first phase transition, essentially, to a first order 
one. The right-hand panel of Fig. [3] refers to this case. 
Figure 2] shows the phase diagram in the region of small 
a (left-hand panel) and over a larger range of separation 
(right-hand panel). The dots connected by the straight 
lines refer to the values of a and T used in Fig. [3] Fi- 
nally, in the left-hand panel of Fig. [5] wc illustrate how 
the condensate changes when both T and a vary linearly 
as in Fig. [3l 



b u 
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Figure 5. Left: Behavior of the condensate when temperature 
and separation change as in Fig. [3] The curves refer to the 
following choices of the parameters: (0) u = p = 1, q — 1/v — 
30; (1) u = p = 1, q = 1/v = 1; (2) u =p = 1, q = l/v = 0.7; 
(3) u = p = 1, q — 1/v — 0.5. Right: Casimir pressure 
(P c = F c /S) for three illustrative cases. Temperature and 
separation change as in Fig. [3] with (1) u = p = 1, q — 1/v — 
30; (2) u = p = 1, q = 1/v = 1; (3) u =p = 1, q = 1/v = 0.7. 
The orange continuous curve refers to the massless case. 

The effects of chiral symmetry breaking on the Casimir 
effect should now be clear. In essence, for a system 
of strongly interacting fermions, it is chiral symmetry 
that controls the fermion mass through the appearance 
of a nonvanishing condensate and, in turn, the parameter 
£ = ma, related to the suppression factor in the Casimir 
force. With the previous results in hands, we have com- 
puted the Casimir force between the layers, and results 
are shown in the right-hand panel of Fig. [5j The thick 
orange curve represents the pressure in the massless case, 
while dashed curves refer to the interacting case with T 
and a varied as before. Because of the temperature de- 
pendence of the condensate, we, in fact, expect two phase 
transitions: one at smaller and one at larger values of a. 
However, effects due to the transition at the smaller value 
of a arc negligible as it can be understood by estimating 
the parameter £ that is bounded by the separation times 
the value attained by the condensate in the infinite vol- 
ume limit, i.e. £ < O(10 ). On the other hand, the 
transition occurring at the larger values of the separa- 
tion leads to sizable effects. Notice, also, that when the 
initial separation is larger (see curve 1 in Fig. [5] (right)), 
suppression of the Casimir force occurs over all the sep- 



aration range. For the case of initial smaller separation 
(see curve 3 in Fig. [5] (right)), the suppression occurs 
over a smaller separation range and the phase transition 
is clearly visible. 

We would like to note here that a behavior similar to 
Fig. has been predicted to occur in bilayer graphene 
with interlayer pairing of electrons and it is expected 
to induce an exotic, possibly experimentally observable, 
form of interlayer superconductivity (3ll | . Whether an ac- 
tual Casimir effect experiment using bilayer graphene (or 
other strongly coupled fcrmionic materials) can be per- 
formed leading to any observable effect of fermion quan- 
tum vacuum fluctuations is currently under investigation. 
Certainly, it seems interesting to examine whether quan- 
tum vacuum energy effects, indirectly, may affect the 
properties of these systems (e.g., stability, separation- 
temperature correlations, etc.), manifesting yet a new 
facet of symmetry breaking, and, eventually, indicate new 
directions of study for the Casimir effect. 

Our goal was to look at the interplay between chiral 
symmetry breaking and boundary effects in the tractable 
and nontrivial case of two parallel layers, for a system of 
strongly coupled fermions. A number of intriguing fea- 
tures arose. First, finite size effects, for the present setup 
tend to change the order of the phase transitions from 
second order (in the infinite volume limit) to first order, 
rendering effects of finite size for the present setup more 
similar to those of a chemical potential than of temper- 
ature. Whether and how an appropriate choice of ge- 
ometry, topology, and boundary conditions may preserve 
chiral symmetry is certainly an important question that 
could be relevant in the context of finite temperature or 
density QCD. Second, we have shown that simultane- 
ous changes in temperature and separation may induce 
an interesting behavior in the thermal dependence of the 
potential and of the condensate that resembles what has 
been discussed earlier for the case of bilayer graphene. 
This similarity provides, in our opinion, motivation to 
consider these issues further. Finally, we have shown how 
chiral symmetry directly influences the Casimir effect, in- 
ducing different phases in the force separated by critical 
lines providing a new example linking the phenomena of 
symmetry breaking with quantum vacuum energy effects. 
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